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FOREWORD 


This report was prepared by Lockheed Missiles and Space Compainy for the 
National Aeronautics and Space Administration under Contract RAS8-26352. 
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the U. S. Army Missile Command. The Contracting Officer's representatives 
were L. A. Kiefling of the MSFC Aero-Astrodynamics Laboratory and J. C. 
Robinson of the Structures Division at Langley. 

The demonstration problems presented in this report were prepared and exe- 
cuted by C. L. Yen and R. A. Moore under the direction of W. D. Whetstone. 
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INTRODUCTION 


A series of examples are presented to indicate some of the principal functions 
of the SPAR system (Reference 1), and to illustrate SPAR’s control card-data 
card structure. For each of the eleven examples, information in the following 
categories is given: 

o A description of the problem and, in most cases, 
comparisons with analytical solutions. 

o A list of the input cards . 

® A printout of the Table of Contents of the 
direct access library into which all SPAR 
output was directed. 

o A few representative plots. 

Many comments are embedded in the input cards. All text to the right of a 
$ symbol is ignored by SPAE’s input decoder. In Problem 1 the input to 
processors TAB and ELD is extensively annotated. Problem 1 also illustrates 
the most general form of applied loading. In Problem 10 the INC and MOD 
commands are extensively used in ELD input. Problem 11 illustrates the use 
of a spectral shift in solving an eigenproblem. 
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1. RING SECTOR 


Solutions were obtained for a segment of a ring in order to explain in de- 
tail various data for processors TAB and ELD, and to demonstrate the use 
of combined loadings in one load case definition in the Q processor. A 6- 
degree sector of a ring was divided into an 8 x 1 quadrilateral finite 
element mesh. The geometrical and material properties of the ring are shown 
on the figure below: 
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Fixed Edge 


Linear and quadratic thermal loads were applied to the ring, as well as a 
combination loading of nodal forces, thermal loads, specified nodal motion, 
surface pressure, and an inertia load. 
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RING sector 


1.5 


2 . VIBRATION OF A CIRCULAR MEMBRANE 


Natural frequencies were obtained for lateral vibration of a simply supported 
circular membrane in uniform tension. A quarter of the membrane was divided 
into 9 X 10 finite elements to obtain the numerical solution. The gecmetri- 
cal and material properties of the membrane are shown on the figure below: 



Numerical results are listed below: 
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3. CIRCULAE PLATE 


The quarter- circle finite element mesh shown on Figure 3-1 was used to solve 
the problems described in Sections 3.1, 3.2, and 3.3. 


B 



F i-S- 3-1 Circular Plate 

3.1 Vibrational Characteristics . With edge AB fixed, and edge OA a symmetry 
plane, solutions were obtained for (l) OB a symmetry plane, and (2) OB an anti- 
symmetry plane. SPAR results are compared with Reference 2, pp UU9, 450. 



OB = 

Symmetry Plane 

OB = 

Anti -Symmetry Plane 


2 

2 


2 


‘“lO 

“^12 

1 

0 

OJ 

3 

■ "11 

SPAR 

.28991 X 10^ 
\ 

.33770 X 10^ 

.43923 X ICp 

.12547 X 10^ 

Ref. 

2 .28638 X 10^ 

.33470 10 ^ 

.43474 X 10^ 

.12371 X 10 ^ 

3.2 

Buckling due to Compressive Edge 

Load. Both OA and 

OB are synmietry 


planes. Unit uniform compression is developed by applying = 1 along free 


3.1 



edge AB. With edge AB then clamped, (critical) is computed. The result 
indicated below is compared with Reference 3, pp 389-392, 

(Critical) 

SPAR .40298 X 10^ 

Ref. 3 .40326 X 10^ 

3*3 Lateral Deflection, iiacluding Prestress Effects . OA and OB are symmetry 
planes. Uniform tension is developed by applying = 100 along free edge 
AB. AB is then clamped, and unit uniform lateral pressure is applied. The 
result indicated below is compared with Reference k, pp 39I-393. 

Lateral Deflection 
at Center 

SPAR .75739 X 10"^ 

Hef. 4 .75960 X 10"^ 
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CIRCULAR Plate 
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buckling mode. critical load = . 4 02798X 1 0"^ 
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4. RECTAKGULAR PLATE 


The 5 by 8 rectangular mesh shewn on Figure U- 1 was used to solve the 
problems described in Sections 4.1-U.9. 


y 



Fig. 4-1 Rectangular Plate 


4*1 Vibra tional Characteristics . Edges AB and BC are simply supported; 
and OA and OC are in symmetry planes. Results are compared with Reference 2, 
p 443, Eq. (188) 


.12906 X 10 


2 2 
"”13 “^31 

.19435 X 10^ .U8406 X 10^ 


2 .12906 X 10^ .19428 X 10^ .48336 X 10^ 


4.1 



Buckling dua to Unidirectional Compression . OA and OC are in symmetry 
planes. With AB free, a unit compressive edge load, N = -1, is applied 
along BC. AB and BC are then simply supported, and the critical value of 
Ny conrputed. Results are compared with Reference 3, pp 351-353 

N (Critical) 

V 

SPAR .52306 X 10^ 

N. 

Bef. 3 .5230I1 X 10^ 

^”3 Effect of Prestress on Buckling Load . OA and OC are in symmetry planes. 
With edge BC free, a tensile load, = 2000, is applied uniformly along AB. 
Effects on lateral stiffness of this constant load are included in the sub- 
sequent analysis. A unit edge load, N^. = -1, is applied along BC. AB and 
BC are then siriply supported and (critical) is computed. Results are 
compared with Reference 3j pp 356-358. 

spar .1270U X lo"* 

Bef. 3 .12700 X 10^ 

.io Uniform C o mpression . OA and OC are in symmetry planes. 
Uniform ccmpression, N . = -1, is applied along AB and BC. AB and 

BC are then simply supported and N (critical) is computed. Results are com- 
pared with Reference 3 > PP 356-358. 

R (Critical) 

■sbar .18830 X iq3 

Bef. 3 .18829 X 1 q 3 

of Pr estres s on Vibrational Characteristics . OA and OC are in 
symmetry planes . 'The plate is uniformly pre- tensioned by applying N = W = 
1000, along AB and BC. With AB and BC simply stj)pbrted, vibrational inodes 



were then computed. 


SPAR 
Ref. 2 


Results are compared with Reference 2, p 435, 443. 


.8x445 X 10^ 
.8x444 X XO*^ 



.46029 X 10^ 
.46022 X XO^ 



.90354 X XO^ 
.90282 X XO^ 


Buckli ng due to Linearly Varying Edge Load . OC is in a symmetry plane. 
With horizontal edge QA and BC free, edge, AB is loaded with N varying 
linearly from N = 1 at A to W = -1 at B. Edges OA, AB, and BC are then 
simply supported and N (critical) computed. Results are compared with 
Reference 3, pp 373-379. . 


SPAR 
Ref. 3 


N (Critical) 

.69434 X 10 ^ 
.653^ X 10 ^ 


Buckling due to Edge Shear . A constant unit shear load, N , is applied 
to the edges. All edges are then simply supported, and the critical shear 
load computed. Results are compared with Reference 3, pp 379-393. 

(Critical) 

SPAR .36190 X 10^ 

I^ef. 3 .36152 X 10^ 


^ckling due to Thermal Load . OA and OC are in symmetry planes. Edges 
AB and BC are simply supported. The critical temperature, T (critical), 
producing lateral buckling is computed. Results are compared with Reference 
pp 356-360, and Reference 10, Chapter l4. 

T (Critical) 

SPAR .43937 X 10 

3 .43935 X 10 


4.3 



OA and OC 


4.9 Effect of Thermal Prestress on Vibrational Characteristics 
axe in symmetry planes. AB and BC are simply supported. Effect 
tlaermal load, T =2, are included in computing vibrational modej, 
are compared with Reference 2, p 435 > P 443. 


s of a 
. Results 


SPAR 
Ref. 2 


.70315 X 10^ 
.70307 X 10^ 



.17156 X 10^ 

.17149 X 10 ^ 
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31 

,44810 X 10 


.44741 X 10 
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5 - VIBRATION OF A COMICAL SHELL 

Vibrational modes of a simply supported truncated conical shell were obtained. 
One quarter of the shell was used for the finite element model. The quarter 
shell was divided into a 9 x 8 finite element network. The seven lowest 
frequencies were obtained. These seven frequencies correspond to half waves 
in the meridional direction with wave numbers three to nine in the circum- 
ferential direction. The geometrical and material properties of the shell 
are indicated in the figure below: 
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Results are compared with Reference 11; 
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Since Reference 11 neglected longitudinal inertia terms, the resulting fre- 
quency values are higher than the exact solution. For low values of circum- 
ferential wave ninriber, most of the ener®r of the shell is associated with 
membrane action. Therefore, neglecting the longitudinal inertia terms re- 
sults in significant error. 
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6. CYLINDRICAL SHELL 

The 10 by 10 cylindrical mesh shown on Figure 6-1 was used to solve the 
problems discussed in Sections 6.1 throvigh 6.5. 



Fig. 6-1 Cylindrical Shell 


6.1 



6.1 Vibrational Characteristics . With edge AB simply supported, vibrational 
modes were computed with (l) BC in a symmetry plane, and (2) BC in an anti- 
symmetry plane. Results are corapsired with Reference 6, p 2154. 


C ir cumf eren tial 
Wave No . , n 


SPAR 


2 


0 ) 


on 


Reference 6 


4 

5 

6 

7 

8 

9 

10 

11 


.1457 X 10 ^ 

.8870 X 10^ 
.6393 X 10^ 

.5828 X 10' 
.6608 X 10c 
.8548 X 10? 
.1168 X 10? 
.1609 X 10° 


.1425 X 10^ 
.8863 X 10^ 
.6330 X 10' 
.5813 X 10' 
.6596 X 10' 

.8493 X 10? 
.1125 X 10? 
.1587 X 10° 


6.2 Buckling due to Axial Compression . A unit compressive edge load, N , 

z 

is applied to AB. AB is then simply supported and critical value of N 

z 

computed. Results are compared with Reference 3, pp 457-458. 


SPAR 
Ref. 3 


N (Critical) 
z 

.6068 X 10^ 

.6075 X 10^ 


6*3 Effects of Pretensioning on Buckling due to External Pressure . The 
shell is pretensioned by a line load, N = 2200, along AB* AB is then simply 
supported, and the critical external pressure, P^, is ccanputed. Results are 
compared with Reference 5, p 425, Figure 9 . 

P 

c 

spar .1415 X 10^ 

PeP' 5 .l406 X 10 ^ 

Buckling due to External Pressure . Edge AB is simply supported. The 
critical pressure, P^, is computed. Results are ccanpared with Reference 5 , 
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p 433, Figure 15. 


P 

c 

SPAR .1084 X 10^ 

Ref. 5 .1045 X 10^ 

^•5 of Pressurization on Buckling due to Axial Compression . With 

edge AB free, an internal pressure, P = 220 , is applied. The stiffening 
effect of this constant preload is included in the subsequent analysis. A 
unit compressive load, is applied along AB. AB is then simply supported, 
and (critical) computed. Results are compared with Reference 4, p 425, 
Figure 9« 

N (Critical ) 

SPAR .6082 X 10^ 

Ref- 4 .6150 X 10^ 



6.3 



#XQT TiB 

.START _ . — 

TITUF* CYLtNBRlCAL SHELL PWOBU^HS 
TEXT* 


GENFWATF basic TABU’S PfFIHjNG STRUCTURE 


THE following PRORLEHS IN CYLINDRICAL ShfUS ARt SOLVfcR IN THIS RUNi 

1 VIBRATION OF A CYLTNORTCAL SHELL 

2 BUCKLING OF A t Yt. I NOR Ic aL SMtLl OUF lO UWIFORH AVIAL COMPRESSION 
buckling of a axially Ppr-SfRESSEO CYLINOLR OUE TO PRESSURE. LOAD, 

4 RUCKLING OF A CVLINORICAl SHELL DUE TO PRESSURF LOAD 

5 buckling of a pressurized CYLINOFR ouF to AXIAI COMPRESSION 

MATERIAL CONSTANTS ^ 

1 , 1+8 .3 .1 . 1 - 4 * 

J D 1 NT.. LOCATION* 

FORMAT= 2 * 

1 10,0 0,0 0 . 10,0 90,0 

U 10,0 0,0 5 * 10.0 90,0 
JOINT REFERENCE FRAME ASSIGNMENT* 

NRFFs-U l» 121 * 

.S.FlKLU_Si:CUaN._PROPER.T.IES .5 

j 0 * 1 s 

constraint CASE IS 

„ S SIMPLY SLIPF'ORTEO AT OWE £NO» SYMMETRY A? 

$ the other end and at planes TKETAsO and 

90 degrees, 


0, 11 IMl 

5 . 


$ 

symmetry PLANt = 3 

symmetry PLAnEsZ 
SYMMETRY' PLANtsl 
ZERO 1,?1 lllil?!* 

constraint case 2* 




symmetry PLANEs? 

. . symmetry plane =2 
antisymmetry PLANfsi 
ZERO l,?t nl,12is 
. C QNS..T R AJ Nl_c. A 8 E _3 i 


SIMPLY ruppohtfo at one end* symmetry At 
J.HE.. other end .and at PLANE...THA.T.A.=.Q.,_A.NT1* 
symmetry at THfcTAsVO DEGREES, 


S 

symmetry PLANEs 3 

SYMMETRY PLANE =2 
symmetry PLANEsl 

-fX(l.t..£LD 

E43 s 

1 Z 13 12 1 10 10* 

. PXQT TDPO .. 

3 

*xoT e 

$ 

«X0T EKS 
S 

PXQT- K 

$ 

9XOT INV 

...S 

PXQT M 

•XOT EIG 
$ 

*• 


symmetry at planes at THETA =0 and 90 
DEGREES, and AT ONE END OF THE CYLINDER, 


.,_REaD....ELEME.NI ..DEFl.NI.TIONS 

, ' 

, analyze element IMTERCONNECTIVITY 

. form element data packets 


.INSERT X* S INTO element DATA PACKETS 
, FORM system K 

, factor matrix in k spar format 

”, FORM system M 

, solve SYSTEM EIGENPROPLEM 

FIND EPFOIJENCIES AND CORRESPONDING MODE 
shapes for free vibration UF the CYLlNOr.R 


itEPRODUClBILITY 


Oi^' x'rAil/ 


L FAGS 2 S FCCF. 


6.U 



•XOT 

•XOT 

s 

% 




...^XQT 

»XOT 

(•Xf)T 

flXOT 

l»XfJT 

„S 

s 

s 

«XOT 

flXDT 

<»XOT 

OXOT 

»XQT 

$ 

S 

VXOT 

fxin 



<»X0T 
$ . 
$ 

J»XQT . 

flXQT 

«KOT 

<»XOT 


UNofR CONSTHAI^T CASt 1 

HESKT IMIT=10,1 ■ 

PWIMT f) 0 0 0 IS 
I NV 

RfcSET C0Ns2$ 

Eir. 

find f RE(Jl(feNCIf..S AND r.ORNFS^’nMDTNfi MOpf 

- shapes rOK FHFt vibration UF THf r.YLlHOr 

UNOF'R cnNSTRAINT CASE 2 

RF.SFT INITttlO.COMri-S 
PRINT 0 0 0 0 JS 

■' » OEFInf APRLUn LOADS 

CASF Jt UNIFOrH axJAL compressive FOftCE apply at one fND 

NODAL . PORCESi MOi'FNTri . 

Ill 3 *>0,7Ssa 2 tos 

ll?3-),b7oflRlS 
CASE ?t untforh radial pressure . ' 
nodal pressure 
1 -1«0 1?1 1$ 

. INV 

RESET CONs3S ~ 

, SOI VE I OAO CASE 1 

RESET L1=1*LHs1»CON=J* 

• form System KOt axial compression 

c I u 

.-.L-Vr;^ - FjNO .HtirKI.ING AXIAL LOAD. DE THE CYLINDER 

RESET PR0BI.EMsSTaB1i.,ITY*IMITs^I»(JPPERs 1 tNDYN=I0.5 

• form linear COMRINATION of MATRTCFS 
FORM stiffness matrix of an axjally pre- 
STRlSSFD CYLINDRICAI, SHELLt AXUL TENSJl I 

, FORCE = 2POO, 

.KtSE.T..RcKXKG»Dlcic.GP = KG«C2 = «2?.00,$ 

INV ■ ~ - •.••• 

reset KKKXKGtCClNa3$ 

DSOl , _ 

RESET L,ls2tL2c2*cON = 3,K=KXKO)i; 

KG ‘ 

RESET . WSE0s2S 

INV 

RESET K=KXKG* 

EIG 

FIND HlicKUING PREISSURF OF PRE-5 TRp SStO 
* ,* CYLiNnEP 

DSo[ ^ ' upper- 1 . PROBsST.AB» NDYNs.1 0» K=1cXKGS 

reset Ll=Pfl.?s2»C0Ns3$ 

RESET OSEOrpS ' * UNIFORM PRESSURE 

EIG 

-- -FJND..R|irKLING. PPE.SSl.iRE,..OF THL CYLINDFR 

reset INITs«.UPPeK= 1 ,PROesSTAB.NOYN=iOS u Y <- 1 N[;r,K . 

LCM 

; ' - ■ form Stiffness matrix of a pressurized 

RFSft D.Kvi.r n, . n, , CVlINdfr* INTERNAL PRFSSURL = 2?0 

RESET R=KYKGfnisK«O?=KG»C2s-?20*S 

,.1NV .. .. ,. ' _ 

RESET KcKYKGf C0N-3S ~ '' " 

DSOI, 

KFSET K = KYKGif.0Na3*( 1 = 1»U251 s 

KK ■ -- ■■■ 

INV 


6.5 



KtSET KsKYKKS 
<»XOT f.IG 

S FIND P.nr.KLI'U; axial LOAP of the P(»r« 

s PMFS.’jHrTZFO r.VLlNDLR 

RFSFT lNlT = «fllPFFRs1 ♦•’ROBsSt ABt NDYtt= 1 0 t KsK VKf.T 
»X0T OCU • FXECUTF UATA roHf'LfX inlunY program 

TOC ts 



TAUle; lip CoNTtNTS* OAI, 1 

cylinuhical Shell ^’RoeLt^-is 


seq 

RR 

date 

TIMfc 

R 

hORoS 

'* 

11 

061174 

121857 

0 

18 

2 

-12 

061174 

121857 

c 

121 

3 

17 

06u74 

121857 

6 

l2 

4 

16 

061174 

12165? 

0 

J9 

5 

19 

061174 

121856 

ft 

14C 

6 

24 

061174 

121658 

c 

iO 

7 

25 

061174 

121901 

0 


a 

38 

06 i i7a 

121901 

c 

121 

9 

43 

061174 

121901 

a 

19 

1 j 

44 

061174 

121902 

c 

121 

n 

49 

06 i 174 

1?19.)3 

ft 

l2i 

12 

54 

06 i i74 

l?l9o5 

<2 

121 

1 3 

59 

09 11 74 

1219o6 

A 

1>89 

1.4 

96 

061174 

121909 

0 

1 6 jO 

IS 

192 

061174 

121909 

0 

tl 

16 

193 

061174 

121910 

0 

l‘j 

17 

194 

061174 

121910 

0 

12 

JS 

165 

001174 

121910 

h 

5 

1 <) 

166 

091174 

l2J‘?iO 

•T 

ft 

0 

2 a 

167 

061174 

121910 

ft 

1 

2 ^ 

177 

0911.74 

1219J0 

0 

1 

22 

176 

091174 

121910 

ft 

1 

25 

179 

061174 

121910 

ft 

1 

24 

lea 

061 i7a 

121910 

ft 

1 

25 

161 

061174 

1219)1 

a 

1 

26 

182 

061174 

121922 

0 

4038 

27 

52o 

061174 

121925 

ft 

12096 

26 

758 

06u74 

12195« 

0 

3()6o0 

29 

1856 

06u74 

lr»l936 

0 

20 

5a 

1659 

061 17« 

122004 

ft 

22400 

3i 

2659 

061174 

122108 

ft 

533i2 

32 

4563 

06il74 

I22i'26 

ft 

224oa 

53 

5363 

061174 

122133 

0 

10 

34 

5384 

061174 

122133 

ft 

72b0 

35 

5624 

06ii74 

122538 

0 

53ii2 

36 

7526 

06 i 174 

122543 

0 

10 

37 

7529 

061174 

122543 

0 

7c^60 

36 

7789 

061174 

123090 

0 

T 20 

39 

7815 

061174 

123040 

0 

l5 

4 a 

7816 

061174 

123^41 

0 

7^6 

<n 

7851 

061174 

123092 

0 

15 

42 

7852 

06i 174 

125092 

0 

0 

45 


09 11 74 

123209 

ft 

1 ^ 

4 4 

-9756 

061174 

123250 

c 

lit, 

45 

-9782 

091174 

123312 

6 

224oO 

46 

-ias62 

091174 

123325 

0 

4 

47 

-1^563 

091174 

123325 

0 

2Vg4 

46 

1 5687 

061174 

125712 

0 

224qO 

49 

11487 

09u74 

123848 

0 

53ii2 

■50 

’13391 

09ii74 

123900 

0 

^2^ 


-13417 

061174 

123916 

0 

22400 

52 

14217 

061174 

124009 

0 

53312 


P0^5 

8LK 

T 

DATA 5lT 


/6LK 

SIZE 

V 


N2 


Ny 

1 

1ft 

0 

JOFI 

6TA0 

1 

e 

l2l 

121 

0 

JREF 

Stab 

2 

6 

1 

12 

1 

aLtr 

Stab 

2 

a 

1 

19 

4 

nOAL 


0 

0 

7 

140 

3 

TfcXT 

6tab 

2 

1 

1 

10 

1 

matc 

6 TAB 

2 


121 

3b3 

1 

JLOC 

6TA6 

2 

s 

121 

121 

ft 

j«rF 

BtA6 

2 

6 

1 

19 

I 

3* 

6TA8 

2 

li 


121 

0 

CUfj 


1 

0 

121 

121 

0 

CUH 


c' 

0 

121 

121 

0 

CON 


i 

0 

1^1 

104 9 

1 

qJJT 

6TAH 

<r 

19 

IOC 

696 

A 

V 

def 

t43 

11 


1 

2 

0 

GO 

E 43 

i 1 

i\ 

1 

15 

jK 

GTir 

E 43 

i 1 

4 

12 

12 

0 

NtLZ 

91 AB 

i 

IJ 

i 

5 

0 

kE 


0 

0 

1 

0 

0 

n5 


0 

(. 

1 

1 

3 

F.L7S 

'''AHL 

ft 

0 

1 

i 

0 

EL7S 

LtVF 

0 

0 

1 

1 

0 

f LTS 

HUOy 

0 

0 

i 

1 

V 

elts 

^SCt 

a 

0 

1 

1 

,ft 

elts 

Ntl-S 

0 

0 

1 

1 

0 

elts 

Lt3 

u 

0 

ii;i 

1 544 

0 

kMap 


541 

29 

l2l 

1344 

0 

a^ap 


1441 


100 

3oft 

a 

E43 


n 

4 

24 

2ft 

0 

OIP 

L 43 

u 

4 

121 

224ft 

1 

K 

Spar 

36 

541 

I2l 

31 3b 

• 

1 

INV 

rt 

1 

l«'4i 

1*21 

2240 

1 

H 

SrAr 

56 

S4i 

10 

Ift 

-i 

viPR 

EvAL 

1 

0 

l2l 

720 

-1 

VIBR 

6 

1 

0 

i'il 

3156 

1 

INV 

N 

2 

I44i 

1C 

10 

-1 

ViBR 

t v AL 

2 

0 

121 

720 

-1 

VlBR 

u 

2 

0 

1*?] 

720 

-1 

c.iFM 


ft 

i 

I 

15 

4 

CASE 


0 

1 

121 

726 

-1 

F H 


u 

2 

1 

15 

4 

CASE 


u 

2 

0 

ft 

0 

LOIR 


1 

2 

1^1 

3156 

1 

INV 

K 

3 

1441 

l2i 

726 

j 

SSUL 

0 

1 

1 


2240 

1 

K6 

SrAr 

36 

54J 

4 

4 

-1 

. 3 TA 0 

EvAl 

1 

0 

121 

726 

-i 

sTab 

6 

1 

0 

l2| 

2240 

1 

kXkg 

SpAR 

36 

54 1 


3156 

1 

IfNV 

Kxng 

3 

144l 

12 1 

720 

1 

SSOL 

U 

1 

2 

l2l 

121 

2240 

3130 

1 

I 

kg 

INV 

SPAH 

NXKG 

36 

1 

54 1 
I44i 


6.7 



53 

•I_6l2l„ 

061174 

12_4J)21 

0 - 

4 

54 

-16122 

061174 

124021 

0 

2Vo4 

55 

16226 

061174 

124306 

0 

7 26 

56 

-16252 

061174 

124323 

0 

224oO 

57 

•17o52 

06ii74 

12432? 

0 

4 

56 

-17053 

061*74 

124327 

0 

29o4 

59 

„ 17157. 

061174 

.1245^7. 

..A_. 


60 

17957 061174 

124602 

0 

533i2 

6l 

1987o 

061174 

124612 

0 

726 


19896 

061174 

124659 

0 

22400 

03 

20696 

061174 

I247j7 

0 

S33i2 


22600 

06ll7a 

124731 

0 

4 

05 

22601 

.06.1 JL74 .12.9731 

0 . 

29fl4 


J21 

726 

121 

4 4 ' 

I 2 l 726 < 

-12J iza). 

12| 3136 

121 726 

I 2 J 2240 

l2l 3136 

4 4 . 

-I2j 726 


i STaB tVAL 
I STA0 

1 sSOL y 

I kG S^'Ar 

a sTas £val 
I SVA9 y 
A KfKG SPA f< 
I r'»V ^KG 
I SSOL ^ 

1 kg Span 
1 jNV KvRG 
I sTA8 £VAL 
i slAiLy. 


i 0 
1 0 
I 2 

it* bH\ 
1 0 
1 0 
3<> b'4i 
5 ia^i 

i 1 

36 541 

1 l« 4 i 
1 0 
1 0 


6.8 





SEQ 


VIBRATIONAL MODE. FREQ C HZ 3 


. 4 024 1 5X 1 0"^ 



CYLINDRICAL SHELL PROBLEM 


SCALE 




6.10 


SEQ 1 


VIBRATIONAL MODE 


FREO CHZD 


. 384223X 1 0 


02 



6.11 


SEO 1 


buckling mode. 


CRITICAL 


LOAD 


- 606797X 1 0 


04 



CYLINDRICAL SHELL PROBLEM p 

scale 


6.12 


SEO 


] 




6.13 


SEQ 


buckling mode, critical load 


.608 1 79X 1 0"^ 



a 


CYLINDRICAL SHELL PROBLEM 


SCALE 




7. BUCKLING OF A CYLINDRICAL SHELL DOE TO TORSIONAL LOADING 

A simply supported cylindrical shell is 'subjected to shear load at the ends 
The shell is divided into ten elements along the axial direction and twelve 
elements along the circumferential direction. The geometrical and material 
properties of the shell are shewn on the- figure below: 


z 



E = .3 X 10® Ibs/in^ 
V = .3 

t = ,1 in 


Results are compared with Ref. 5, p h48. Figure 21. 

N^e (Critical) 

SPAR .3437 X 10^ 

Ref. 5 .3423 X 10^ 
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SEQ 1 


Buckling mode, critical l □ a d = . 343730x1 



Shear buckl I ng of a cylindrical shell p 

' scale 
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8. SPHERICAL SHELL 

The 6 by 10 mesh shown on Figure 8-1 was used to solve the problems dis 
cussed in Section 8.1 through 8.5 




Section 

a degrees 

t in . 

E Ib/in^ 

V 

8.1 

8.2 

60 ° 

5. 

5. 

10 !^ 

• 3 

8.3 

8.4 

35 

10° 

3. 

.396 

10 ' 

♦ J 

1/6 

1/3 


Fig. 8-1 Spherical Shell 


Mas s ' o li 

Density Ib-secVin^ 

.1 


.1 
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8.1 Vibrational Characteristics of a 60° Sector . With edge AB simply sup- 
ported, vibrational mode? are computed. Results are compared with Reference 6, 
p 1363 , Table 2. 


SPAR 
Ref. 6 


U), 


01 


.9209 X 10 
.9197 X 10^ 


.1736 X 10 ^ 
.1764 X 10 ^ 


8.2 Vibrational Characteristics of a 30° Sector . With edge AB fixed, 
vibrational inodes are computed. Results are compared with Reference 6, 
P 1363 . 


SPAR 
Ref. 6 


.4743 X 10^ 
.4623 X 10 ^ 


8-3 Stresses due to Pressure Loading . With edge AB fixed, unit internal 
pressure was applied. Results are compared with Reference 4, pp 553-554, 
Figure 273-274. 


SPAR 
Ref. 4 


Hoop Stress Meridional 

at Apex Moment at AB 


48.0 35.8 

47.5 36.6 


Buckling of 10° Sector due to Pressure . To obtain a state of uniform 
compression, unit pressure, P, was applied with displacements -hormar” to > 
conical surface ABC constrained along AB. Edge AB was then fixed, and 

P (critical) computed. Results are compared with Reference 8, p 450, 
Figure 6. 

P (Critical) 

spar .236 X 10 ^ 

7 .234 X 10^ 
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7 

7 
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7 1 
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.7 
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$ 

S 
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1 
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SIMPLY SUPPORTED AT EDGES* SYMMETRY AT 
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.2._9.-10...1 l_9_6S 

E33 $ 

1 2 3 2 6$ 

..... 9XOT TOPO - 

$ 

• XQT E 

$ 

^XOT EKS 
S 

."PXQT M 
$ 


• ANALY7E ELEMENT INTERCONNECTIVITY 
, form elfhfnt data packets 

• insert K* S into element data PACKETS 

• form SysIEM h 
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»XQT INV 
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$ 
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S 

RESEI_I.NJ.tsBS 
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«XOT 
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. form system K 
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JOINT LOCATIONS, 30 DEGREE DOME 
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0,1679 

0,670? 

1 ,5071 
- 2,6734._. 
4,1655 
5,9770 
8,1035 
10,5347 
13, 2624 
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— FIND DTSPLACEHENT -PE the 3'i DEGREE DOME 
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^ FIND STRESSFS OF THE OOMt 

. print stre:ss data 

_L_£P.IW.t_SXK£.SStS. .OF_lHL.bOME . 


JOINT LPCATIONSi 10 OFGREE DOME 


1 

0.0 

0.0 

0,0$ 





2 - 

-.1,7452 n^n 

-0,0152. 

t-,.7.4 5 2-90 ,i3 

. 0-*.0 1 . 52 . 

7_ 

1 £ 

9 

3.4900 
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9- BEAM 

The 10 element mesh shown on Figure 9~1 was used to solve the problems dis- 
cussed in Section 9-1 throTigh ^. 8 . 
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Fig. 9-1 Beam 
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9.1 Vibration of a Cantilever . With point A fixed and the beam constrained 
to move in the yz plane, vibrational modes were computed. Results are ccmr 
pared with Reference 1 , p 339 - 
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9«2 Column Buckling . With point A fixed, the lowest critical value of, 
compressive load at B was computed. 

.98695 X 10^ 

.98695 X 10^ 

9 ‘3 Beam-Column Effects . A compressive pre-load of 1000 lbs was applied. 
With simple support at A and B, lateral displacement under a 10 lb lateral 
load at C was computed . 

SPAR 
Ref. 8 

9 -^ Effects of Pre-Tension on Vibration . A tensile pre-lo^ of 1000 lbs was 
applied. With simple support at A and B, vibrational modes were computed. 
Results are conqpared with Reference 2 , pp 374 - 376 . 
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value of an x-directxon point force at B was computed. Results are compared 
with Reference 3j PP 259-260. 
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9*6 Buckling due to Pure Bending . Critical values of equilibrating point 
moments (about the y axis) applied at A and B are computed. Results are 
compared to Reference 2, pp 253-257. 
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9-7 Buckling due to Thermal (|joad . With A and B simply supported, the criti- 
cal temperature was computed. 
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9.8 Effects of Thermal Pre~Stress on Vibration . A 10° temperature load was 
applied to the beam, with A and B simply sr^sported. Vibrational modes were 
conpated. Results are compared with Reference 2, pp 37^-375 
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10. LUT (SATURIf 5 lAUWCHER UMBILICAL TO^R) 

Inertia prestress vibrational modes were calculated for the LUT. This 
problem illustrates use of the MOD and INC commands in processor ELD. 

The LUT is a multi-story frame structure, the floors of which are supported 
by columns and diagonals. Non- structural mass was attached at the four 

of each floor as rigid lumped mass. Varying cross-sectional proper- 
ties for the general beam sections were defined throughout the structure. 

The accompanying figures and table show the section property data set dis- 
tribution. The geometrical and material properties for the LUT are listed ' 
below: 


E = 3.0 X lo"^ Ib/in^ 


V = 0.2998 

Specific _ oA IV /. 3 

WeigDt ■ 




Fixed 

Supports 


60' 111’ 2" |-ii 

A downward inertia load of 386. 088 (equivalent to a one-g loading) was used 

as a prestress load, and cantilevered vibrational modes were calculated for 
the prestressed structure. 
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Table 10.1 Section Data Sets for Floors and Bays 
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The elastic modes are stunmarized below: 
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UNDEFORMED STRUCTURE 
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11, MECHANISM 


Free-free vibrational modes were obtained for a linkage in order to demon- 
strate the application of a spectral shift in the eigensoliition process. 

The linkage was defined as a syston of two beam segments connected through 
a zero-length 6x6 directly specified stiffness matrix. Planar motion was 
specified for the otherwise unconstrained configuration. The geometrical 
and material properties of the linkage are shown on the figure below: 
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Eigenvalues were desired in the neighborhood of the value c = 0.1 x 10 , which 

was selected as the shift point for the eigenvalue solution of K-cM (see Reference 1, 
Section 3-17) ♦ Four "zero" eigenvalues were obtained, representing the three 
rigid body motions and a mechanism mode. 

The elastic modes are summarized below: 
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11. MECHANISM 


Free-free vibrational modes were obtained for a linkage in order to demon- 
strate the application of a spectral shift in the eigensolution process. 

The linkage was defined as a system of two beam segments connected through 
a zero- length 6x6 directly specified stiffness matrix. Planar motion was 
specified for the otherwise unconstrained configuration. The geometrical 
and material properties of the linkage are shown on the figure below: 
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Eigenvalues were desired in the neighborhood of the value c = 0.1 x 10^^, which 
was selected as the shift point for the eigenvalue solution of K-cM (see Reference 1 
Section 3-17) - Four "zero” eigenvalues were obtained, representing the three 
rigid body motions and a mechanism mode. 
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The elastic modes are summarized below: 
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